Integrative multi-dimensional characterization of striatal projection neuron heterogeneity in adult brain.

Striatal projection neurons (SPNs) are traditionally segregated into two subpopulations expressing dopamine (DA) D1-like or D2-like receptors. However, this dichotomy is challenged by recent evidence. Functional and expression studies raise important questions: do SPNs co-express different DA receptors, and do these differences reflect unique striatal spatial distributions and expression profiles? Using RNAscope in mouse striatum, we report heterogenous SPN subpopulations distributed across dorsal-ventral and rostral-caudal axes. SPN subpopulations co-express multiple DA receptors, including D1 and D2 (D1/2R) and D1 and D3. Our integrative approach using single-nuclei multi-omics analyses provides a simple consensus to describe SPNs across diverse datasets, connecting it to complementary spatial mapping. Combining RNAscope and multi-omics shows D1/2R SPNs further separate into distinct subtypes according to spatial organization and conserved marker genes. Each SPN cell type contributes uniquely to genetic risk for neuropsychiatric diseases. Our results bridge anatomy and transcriptomics to offer new understandings of striatal neuron heterogeneity.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning.

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Molecular rhythm alterations in prefrontal cortex and nucleus accumbens associated with opioid use disorder.

Severe and persistent disruptions to sleep and circadian rhythms are common in people with opioid use disorder (OUD). Preclinical evidence suggests altered molecular rhythms in the brain modulate opioid reward and relapse. However, whether molecular rhythms are disrupted in the brains of people with OUD remained an open question, critical to understanding the role of circadian rhythms in opioid addiction. Using subjects' times of death as a marker of time of day, we investigated transcriptional rhythms in the brains of subjects with OUD compared to unaffected comparison subjects. We discovered rhythmic transcripts in both the dorsolateral prefrontal cortex (DLPFC) and nucleus accumbens (NAc), key brain areas involved in OUD, that were largely distinct between OUD and unaffected subjects. Fewer rhythmic transcripts were identified in DLPFC of subjects with OUD compared to unaffected subjects, whereas in the NAc, nearly double the number of rhythmic transcripts was identified in subjects with OUD. In NAc of subjects with OUD, rhythmic transcripts peaked either in the evening or near sunrise, and were associated with an opioid, dopamine, and GABAergic neurotransmission. Associations with altered neurotransmission in NAc were further supported by co-expression network analysis which identified OUD-specific modules enriched for transcripts involved in dopamine, GABA, and glutamatergic synaptic functions. Additionally, rhythmic transcripts in DLPFC and NAc of subjects with OUD were enriched for genomic loci associated with sleep-related GWAS traits, including sleep duration and insomnia. Collectively, our findings connect transcriptional rhythm changes in opioidergic, dopaminergic, GABAergic signaling in the human brain to sleep-related traits in opioid addiction.

Transcriptional Alterations in Dorsolateral Prefrontal Cortex and Nucleus Accumbens Implicate Neuroinflammation and Synaptic Remodeling in Opioid Use Disorder.

BACKGROUND: Prevalence rates of opioid use disorder (OUD) have increased dramatically, accompanied by a surge of overdose deaths. While opioid dependence has been extensively studied in preclinical models, an understanding of the biological alterations that occur in the brains of people who chronically use opioids and who are diagnosed with OUD remains limited. To address this limitation, RNA sequencing was conducted on the dorsolateral prefrontal cortex and nucleus accumbens, regions heavily implicated in OUD, from postmortem brains in subjects with OUD. METHODS: We performed RNA sequencing on the dorsolateral prefrontal cortex and nucleus accumbens from unaffected comparison subjects (n = 20) and subjects diagnosed with OUD (n = 20). Our transcriptomic analyses identified differentially expressed transcripts and investigated the transcriptional coherence between brain regions using rank-rank hypergeometric orderlap. Weighted gene coexpression analyses identified OUD-specific modules and gene networks. Integrative analyses between differentially expressed transcripts and genome-wide association study datasets using linkage disequilibrium scores assessed the genetic liability of psychiatric-related phenotypes in OUD. RESULTS: Rank-rank hypergeometric overlap analyses revealed extensive overlap in transcripts between the dorsolateral prefrontal cortex and nucleus accumbens in OUD, related to synaptic remodeling and neuroinflammation. Identified transcripts were enriched for factors that control proinflammatory cytokine, chondroitin sulfate, and extracellular matrix signaling. Cell-type deconvolution implicated a role for microglia as a potential driver for opioid-induced neuroplasticity. Linkage disequilibrium score analysis suggested genetic liabilities for risky behavior, attention-deficit/hyperactivity disorder, and depression in subjects with OUD. CONCLUSIONS: Overall, our findings suggest connections between the brain's immune system and opioid dependence in the human brain.

Addiction-Associated Genetic Variants Implicate Brain Cell Type- and Region-Specific Cis-Regulatory Elements in Addiction Neurobiology.

Recent large genome-wide association studies have identified multiple confident risk loci linked to addiction-associated behavioral traits. Most genetic variants linked to addiction-associated traits lie in noncoding regions of the genome, likely disrupting cis-regulatory element (CRE) function. CREs tend to be highly cell type-specific and may contribute to the functional development of the neural circuits underlying addiction. Yet, a systematic approach for predicting the impact of risk variants on the CREs of specific cell populations is lacking. To dissect the cell types and brain regions underlying addiction-associated traits, we applied stratified linkage disequilibrium score regression to compare genome-wide association studies to genomic regions collected from human and mouse assays for open chromatin, which is associated with CRE activity. We found enrichment of addiction-associated variants in putative CREs marked by open chromatin in neuronal (NeuN+) nuclei collected from multiple prefrontal cortical areas and striatal regions known to play major roles in reward and addiction. To further dissect the cell type-specific basis of addiction-associated traits, we also identified enrichments in human orthologs of open chromatin regions of female and male mouse neuronal subtypes: cortical excitatory, D1, D2, and PV. Last, we developed machine learning models to predict mouse cell type-specific open chromatin, enabling us to further categorize human NeuN+ open chromatin regions into cortical excitatory or striatal D1 and D2 neurons and predict the functional impact of addiction-associated genetic variants. Our results suggest that different neuronal subtypes within the reward system play distinct roles in the variety of traits that contribute to addiction.SIGNIFICANCE STATEMENT We combine statistical genetic and machine learning techniques to find that the predisposition to for nicotine, alcohol, and cannabis use behaviors can be partially explained by genetic variants in conserved regulatory elements within specific brain regions and neuronal subtypes of the reward system. Our computational framework can flexibly integrate open chromatin data across species to screen for putative causal variants in a cell type- and tissue-specific manner for numerous complex traits.

Impact of acute lymphoblastic leukemia induction therapy: findings from metabolomics on non-fasted plasma samples from a biorepository.

INTRODUCTION: Acute lymphoblastic leukemia (ALL) is among the most common cancers in children. With improvements in combination chemotherapy regimens, the overall survival has increased to over 90%. However, the current challenge is to mitigate adverse events resulting from the complex therapy. Several chemotherapies intercept cancer metabolism, but little is known about their collective role in altering host metabolism. OBJECTIVES: We profiled the metabolomic changes in plasma of ALL patients initial- and post- induction therapy. METHODS: We exploited a biorepository of non-fasted plasma samples derived from the Dana Farber Cancer Institute ALL Consortium; these samples were obtained from 50 ALL patients initial- and post-induction therapy. Plasma metabolites and complex lipids were analyzed by high resolution tandem mass spectrometry and differential mobility tandem mass spectrometry. Data were analyzed using a covariate-adjusted regression model with multiplicity adjustment. Pathway enrichment analysis and co-expression network analysis were performed to identify unique clusters of molecules. RESULTS: More than 1200 metabolites and complex lipids were identified in the total of global metabolomics and lipidomics platforms. Over 20% of those molecules were significantly altered. In the pathway enrichment analysis, lipids, particularly phosphatidylethanolamines (PEs), were identified. Network analysis indicated that the bioactive fatty acids, docosahexaenoic acid (DHA)-containing (22:6) triacylglycerols (TAGs), were decreased in the post-induction therapy. CONCLUSION: Metabolomic profiling in ALL patients revealed a large number of alterations following induction chemotherapy. In particular, lipid metabolism was substantially altered. The changes in metabolites and complex lipids following induction therapy could provide insight into the adverse events experienced by ALL patients.